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Abstract 

The “classic” analogy of classical repulsive interactions via exchange of particles is revisited 
with a quantitative model and analyzed. This simple model based solely upon the principle of 
momentum conservation yields a nontrivial, conservative approximation at low energies while also 
including a type of “relativistic” regime in which the conservative formulation breaks down. Sim¬ 
ulations are presented which are accessible to undergraduate students at any level in the physics 
curriculum as well as analytic treatments of the various regimes which should be accessible to 
advanced undergraduate physics majors. 
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I. INTRODUCTION 


Countless students in introductory physics learn that the “exchange of virtual parti¬ 
cles” is responsible for the fundamental forces of nature. Several popular introductory 
textbooks contain diagrams which sketch how classical particle exchange could plausibly 
explain the qualitative nature of repulsive forces.^ Furthermore, some texts even attempt 
to construct analogies for how attractive forces could arise from complicated exchanges of 
classical objects.® In this paper, we wish to address the gaping hole in the literature re¬ 
garding how such pictures may be quantitatively useful in understanding the connection 
between fundamental interactions and momentum transfer through mediating particles. 

Just as physical theories are only useful within certain domains of validity, analogies are 
only helpful until their meanings are stretched to a point at which the usefulness breaks 
down. To properly analyze fundamental interactions, the methods of quantum held theory 
provide the tools necessary for obtaining quantitatively accurate results. Ref. [5] provides a 
particularly illuminating discussion of how gravitational, electrostatic and nuclear potentials 
arise as either attractive or repulsive interactions by using the path integral formualtion of 
quantum held theory. Additionally, by casually invoking the energy-time version of the 
Heisenberg uncertainty principle, one may obtain surprisingly accurate information regard¬ 
ing the force laws resulting from electromagnetic and nuclear interactions.® The focus of the 
present work is not to require an idealized analysis within classical mechanics to describe 
the nature of fundamental interactions, but to explore how ehective forces between particles 
which are spatially separated can arise within classical dynamics. 

A student needs only very basic tools to explore the implications of a particular particle 
exchange model. With easily acquired numerical results, an advanced student may apply 
the mathematical analysis required to obtain both exact and asymptotic results. The goal 
of the present work is to present a quantitative approach, accessible at both introductory 
and advanced levels, which thoroughly analyzes a particular model for interactions based on 
classical physics. 

In particular, we consider a system of two massive particles, each of mass M, which 
interact with each other via the exchange of two mediating particles, each of mass m <C 
M, which are taken to always move at speed c and interact with the heavier particles 
through inelastic collisions, always emerging with speed c relative to a stationary lab (or 
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“ground”) frame. Though this model is admittedly artihcial compared to the quantum 
held theories describing the known fundamental interactions, the reasoning required for a 
careful, quantitative analysis are quite useful in understanding the realistic interactions that 
do occur in nature through mediating quantum helds.^ 

A notable shortcoming of the classical particle-exchange analogy is its inability to describe 
attractive forces.^ While it is possible to invoke quantum huctuations in energy to explain 
attractive nuclear forces in a qualitative manner,!^ we emphasize that attractive interactions 
emerge naturally from classical scalar held theory.l^Such a rigorous discussion of the origin of 
attractive forces implicitly requires a discussion of quantum theory, as these interactions rely 
on the wave-like nature of matter. Consequently, such treatment is beyond the scope of the 
present work, as we wish to present a model which may be thoroughly analyzed classically. 

This paper is arranged as follows: in Sec. [IT] we present a model for classical particle 
exchange and explore some basic consequences through simulations and physical reasoning, 
both of which are appropriate for students in introductory physics courses. Sec. |III| contains a 
thorough analysis of the model employing advanced physical reasoning and special functions 
to verify the speculative results obtained through careful estimation in Sec. |TTj Finally, we 
summarize the results in Sec. IIVI 

II. MODEL 

We wish to investigate the classical picture of particle exchange as a mechanism for inter¬ 
action between two massive particles. We imagine two particles each of mass M exchanging 
small particles, each of mass m M as shown in Fig. The analogy is often made to a 
pair of ice skaters (or rollerbladers) tossing a ball back and forth.M^l Each time one skater 
catches the ball and throws it back, a small amount of momentum is imparted to the skater, 
resulting in an effective repulsive force between the skaters which is mediated by the ball 
being tossed. We construct a quantitative model for this type of interaction by taking the 
smaller particle’s velocity to be a constant, given speed c. We choose the label c with no 
reference whatsoever to the speed of light, though we will see that our c plays a role in 
our model which is rather similar to that of the actual speed of light in electromagnetism, 
allowing us to explore a sort of “non-relativistic” limit of the model for speeds n <C c. In 
order to keep the system’s center of mass at rest, we shall consider a symmetric setup in 
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which two small particles are exchanged. When the smaller, mediating particles approach 
each other we assume that they pass through one another without interaction or collide 
elastically.^ 



FIG. 1: Two particles of mass M experience a repulsive “force” which is mediated by the exchange 
of a smaller particle of mass m <C M. 


Since the mediating particles always move at speed c, the collisions involving the massivd^ 
particles with the mediating particles would not result in momentum transfer if the collisions 
were elastic. To obtain nontrivial momentum transfer, we must consider inelastic collisions 
which result in an incremental increase in the system’s kinetic energy after each collision. We 
shall explore whether the work required for this change in kinetic energy may be associated 
with an effective potential energy for the system. Taking the large, right-moving particle to 
be moving at speed v, momentum conservation applied to a single collision gives 

Mvn + me = Mvn+i ~ me, ( 1 ) 

or 6v = Vn+i — Vn = 2^. With repeated collisions of this form, the two massive particles 
will accelerate away from their common center of mass in a manner qualitatively similar 
to the motion experienced by two like charges placed near each other and released. We 
employ two approaches to investigate the quantitative nature of this effective force law. 
First, we simulate the system as described, obtaining numerically an effective force law 
which decreases as for small velocities v e, where r is the instantaneous separation 
between the two massive particles. Second, the discrete sequence of collisions leads to a 
recursion relation which allows us to obtain a closed-form expression for r„, the separation 
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distance immediately preceeding the collision. While exact, this closed-form expression 
for is less than transparent regarding the physics of the system. In the following section, 
we apply continuum approximation to uncover the effective dynamics analytically in various 

limits. 


A. Full simulation 


The full simulation consists of integrating the Newtonian equations of motion for free 
particles moving at constant speeds and monitoring for a “collision” at which point each 
massive particle is given a boost in speed 5v = 2mc/M and the mediating particles are 
reflected with equal momenta in the opposite directions. Letting denote the 

position of the right-moving (left-moving) particle and its velocity, we consider 

the following initial conditions: 

a;d)(o) = -a;(2)(o) = (2) 

yd)(0) = -i;(2)(0) = 0. (3) 


The mediating particles are initially located at the origin and begin moving in opposite 
directions toward the massive particles at t = 0 with speed c. Letting the positions of the 
mediating particles be given by for i = 1, 2, it is an instructive exercise to numerically 
integrate the equations of motion 


dt 

dt 



Vd), 


(4) 

(5) 


with = +c and = —c at t = 0. To monitor for collisions, at each time step At we 
check for the following condition: 


|a.»_xb)|<e, (6) 

indicating that the mediating particle nearest the particle has come within a small 
distance e of the massive particle’s location. When this occurs, we make the following 
adjustment to the equations of motion: 

vii) ^ yd) + ^sign , 

V(t) 
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(7) 

( 8 ) 



indicating that a collision has occurred, resulting in momentum transfer. Results are gen¬ 
erally insensitive to the time-step size, provided e ^ cAt. Fig. depicts the numerically 
computed average acceleration as a function of separation distance for m = 0.005M. For 
the computation of acceleration, we only use the separation distance and corresponding time 
just after collision events, since each massive particle’s acceleration is formally zero between 
collisions. Note that for position measurements which are taken at unequal time increments, 
we require the following discrete represent at iorP^ of its second temporal derivative 


dt'^ 


r=rn 


rn+l-Tn _ Tn-rn-l 
^n + l ^n — 1 




(9) 



FIG. 2; Numerically computed acceleration plotted against separation distance for the right-moving 
mass with logarithmic scales on axes. Also shown is a linear regression for the logarithmic data. 


A strong linear trend on a log-log plot demonstrates the power-law nature of the force 
law. 


dt"^ 


oc r 


hi 


( 10 ) 


with hi ^ —1. This result is consistent with a rough estimation of the rate of momentum 
transfer for u <C c. Each collision is associated with transfer of momentum 


Sp = 2mc. 


( 11 ) 
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For f <C c, the massive particles do not move appreciably during one collision cycle. Let r 
denote the instantaneous separation distance between the massive particles. Beginning with 
the mediating particles at the origin, one cycle requires each particle to cover a distance | 
to collide with the massive particles and then another distance of | to return to the origin. 
Thus, a single collision cycle associated with a momentum transfer 6p requires a time 


6t = 


r 

c 


The average force experienced by each massive particle is then 


TT' _ _ 

^ ave — r-» — 

dt 


6p 2mc^ 


( 12 ) 


(13) 


The validity of this crude estimate will be examined more carefully in the next section, but 
for now it serves to make the results in Fig. appear rather plausible. One might worry 
about the implications of an inverse-linear force law, since this could potentially be associated 
with a logarithmic potential energy function, just as in the case of two uniformly charged 
wires of inhnite length.l^^ In the case of point particles, the potential does not asymptotically 
approach a constant value at large distances and should result in ever-increasing speeds as 
the massive particles move farther away from each other. This does not appear consistent 
with the model, as the mechanism for momentum does not allow the mediating particles to 
travel faster than speed c, so the speeds of the massive particles should be bounded by this 
limit. The resolution of this apparent paradox will be addressed below where we must rehne 
the simulation method in order to access much longer times. 


B. Calculation of collision times 

The results so far suggest a disconnect between the low-energy behavior of the model 
and the high-energy “speed limit” of c, which should be enforced by the mediating particles. 
To obtain some resolution, we must explore extremely large timescales, thus allowing the 
massive particles to approach high speeds, ~ c. Because the time between subsequent 
collisions grows at an accelerated rate as the massive particles spread apart and speed up, 
the basic scheme outlined above becomes impractical. In fact, most of the computation is 
entirely unnecessary since all particles move with constant velocities until a collision occurs. 
Starting from one collision event, the time for the next collision may be computed using the 
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instantaneous velocities of all particles, and this process may be repeated. Though the time 
between collisions grows rapidly, the computation time of this scheme grows linearly with 
number of collisions, not with the elapsed time as before. 

To proceed, let us consider a single collision event shown in Fig. [Tj With both mediat¬ 
ing particles at the origin and instantaneous separation between the outwardly moving 
massive particles, the next collision will occur after the mediating particles have reached the 
massive particles, requiring a time 


5tn 


C-Vn 


(14) 


corresponding to traveling a distance of ^ with speed c — Vn relative to the outwardly 
moving, massive particles. After time Stn has elapsed, collisions occur resulting in the 
mediating particles reversing directions and 


Vn Vn+1 =Vn + 


2mc 

~W' 


(15) 


The cycle completes when the mediating particles return to the origin. By symmetry, this 
also requires time 6tn, so the entire elapsed time for a complete cycle is 26tn, or 

tn+l = tn -\ -• (16) 

C-Vn 

To update the positions of the massive particles, we note that before the collision, each 
particle was moving away with speed Vn with respect to the ground for time 6tn- After 
the collision, each particle moves away from the system’s center of mass for time 6tn with 
the updated speed, Vn+i- Thus, the separation distance increases by an amount 2vnStn + 
2Vn+lStn, or 

f'n+l T 2VnStn T 2Vn-\-lStn- (l^^) 


Eqs. (@-(0 constitute a closed recursion relation which may be iteratively advanced to 
obtain the velocity, separation distance and time corresponding to the beginning of each 
collision cycle. 


For a point of comparison, we may take the approximate force law in Eq. (13) and write 
Newton’s second law for the motion of the right-moving particle, 

2mc? 


M- 


dP 


(18) 


Applying the symmetry of the system, we have r = 2x^^'^ and may change variables, 

d?r Id /.nx 

-= -=- (r^) . 

dP dP 2 dr ^ ^ 

8 


(19) 










FIG. 3: Long-time, large-distance behavior of massive particle speed (blue circles) 
Eqs. (|14[)-([l7l) and compared to low-speed, non-relativistic (NR) approximation in 
provides excellent agreement with the simulation for u <C c. 


computed from 


Eq. (20), which 


Writing r = 2n, where v represents the speed of each massive particle, we may integrate 
both sides to obtain 


v‘^{r) = vl + 


2mc 


M '"7,' 


( 20 ) 


which represents a statement of conservation of energy with a potential energy given by 

f/(r) = ^l„“. (21) 


for some arbitrary length scale a. We refer to Eq. (20) as the non-relativistic approximation, 
as its derivation relies on assuming n -C c. The term “non-relativistic” (NR) as used here 
does not refer to speeds much less than the actual speed of light but those significantly 
smaller than the mediating particle speed c. The role played by c in this model is similar 
to that of the actual speed of light in electrodynamics, but we stress that special relativity 
and the actual speed of light play no role in this model. Improvements to this low-energy 
approximation will be explored in the next section, but we are in a position to compare its 
predictions to the full simulation. Fig. [^depicts the predictions of Eq. (20) compared to the 
actual simulation information contained in Eqs. (p!^-(17). As expected, the non-relativistic 
approximation breaks down as the massive particles’ speeds approach c. For large separation 
distances, the massive particle speeds do not increase as sharply with increasing distance as 
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the non-relativistic approximation predicts. Indeed, once the massive particles reach a speed 
of c, the mediating particles, also traveling at speed c, are unable to catch up to the massive 
particles. Correspondingly the recursion relations break down and no more collisions are 
found. Specihcally, as —>■ c from below, we have 6tn —)■ oo. If the massive particle speed 
becomes exactl}^ c, 6tn does not exist and no further collisions occur. Another possibility 
is that a single collision changes from just below c to just above c. In this case, 5tn 
formally becomes negative and we conclude similarly that no further collisions occur. 

The behavior of the system explored thus far can be summarized as follows: for arbitrary 
initial separations, the massive particles are repelled from each other by the effective force 
provided by the mediating particles. At long times, the speeds (with respect to the ground) of 
the massive particles approach c, the speed of the mediating particles. While an approximate 


statement of energy conservation has been derived (see Eq. (20)) for low speeds u -C c, the 
associated potential is problematic as it has no lower bound for r ^ oo. An unlimited 
amount of potential may be converted into the massive particle’s kinetic energy resulting in 
the erroneous prediction that for any initial separation, both massive particles will continue 
to accelerate rather than asymptotically approach hnite speeds. That the initial separation 
distance has no effect on the hnal speeds of the massive particles suggests that the system 
is not conservative. In the next section, we will carefully examine this system using analytic 
tools to quantitatively explore some of these issues. 


III. ANALYTIC APPROACH 


Exact solution to recursion relation 


The discrete sequence of collisions described by Eqs. (0-(|Irl) can be analyzed exactly. 


yielding a closed-form expression for r„, the separation distance after n collisions. Eq. (15) 
simply states that the velocity increases by a constant amount after each collision, or 

2mnc 


— 


M 


( 22 ) 
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Inserting Eq. (22) into Eq. (17) and using Eq. (14), we have 


Tn+l = rn +2 


1 + 


me 
Vn H-TV 


^ ' 71 /f ’ 

Mi C — Vn. 

2m(n+l) 


M 


2mn 

M 


Tr. 


Proceeding iteratively, 


ri = 


r2 = 


^ 2m\ 

(i + t)(i + ^) 

(i-t) 




Tr,, = 


/ L._n 


M I 1 — ^ 

^ fc=0 \ M 


By employing the Gamma function, which satisheJ^^ 


r(a; + 1) = a;r(a;). 


(23) 

(24) 

(26) 

(26) 

(27) 

(28) 

(29) 


and reduces to the factorial for integer arguments, n\ = r(n + 1), we may write this as 


Tr,. = 


[r(lh]' 


2mn\^ 

1 I ’^o- 


M 


(30) 


The derivation of Eq. (30) from Eq. (28) requires use of Eq. (29), the property,' 


na 


r(x)r(i — x) = 


TT 


and their mathematical offspring. 


r(x)r(-a:) = 


sm vra; 


TT 


xsin vra; 


(31) 


(32) 


B. Limiting cases 


As an exact, closed-form solution, Eq. (30) contains all of the physics we have encountered 


up to this point. The low-energy force law in Eq. (13) was previously derived using physical 


reasoning, but we can demonstrate that it also follows from the exact solution rather than 
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appealing to comparisons such as Fig. To this end, let us define a = ^ and take the 


natural logarithm of Eq. (30), obtaining 


In — = In r (a — n) + In r (a — n) — 2 In r (a) 
ro 

2 " 


In 


-'S' 


(33) 


To investigate the dynamics for m M and n <C c, we examine the limit a —)■ cx) with 
n a. We first apply Stirling’s approximation!^ to the Gamma functions. 


lnr(a±n) ~ (a ± n) In [a ± n], 
lnr(Q;) ~ a In a. 


(34) 

(36) 


Applying the limit n <C a and expanding the logarithms according to 


X 


we recover the result 


(1 ± x) In [1 ± x] ~ a: + —, 


r n 

In — ~ —, 
ro a 


(36) 


(37) 


which is equivalent to Eq. (20) with uq = 0 upon the identification n ^ (see Eq. (22)) 


Alternatively, we may consider the limit n —)■ c. Note that Eq. (30) diverges as n —>■ a. 


indicating that this only occurs as r —)■ oo. Implicit in this relation is the upper limit on 
number of collisions before the massive particles reach terminal velocity, 

M 




2m 


(38) 


We may probe the system at long times by letting n = a — e for e 1. Eq. (30) then 
becomes 

r* F' 9/n/^ Of: 

(39) 


— = (£)■-. 
ro [T (a)] a 


Employing the small-argument expansioiJ^ 


r(e) = --7 + 0(6). 


(40) 


where 7 ~ 0.577 is the Euler-Mascheroni constant, we may expand Eq. (39) to obtain 


r 2 r ( 2 a) 
_ _ 2 __ 

ro a [r(a)]^ 


(1 - 7e). 


(41) 
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Taking e —)■ 0 is equivalent to letting n —)■ c, and we obtain 


4mr (^) 

r -^rr_ = -^ , 0 ^0 


as n —)■ c. 


(42) 


M[rS)]' 

For r > Tc, each massive particle is moving at the same speed as the mediating particles and 
experiences no subsequent collisions with the mediating particles. Some mystery may be 


removed from Eq. (42) by taking the natural logarithm of both sides and applying Stirling’s 


approximation, this time keeping several terms 

1 X 

In r(x) ~ a; In a; — X-In —. 

2 

When the smoke clears, we have the compact result 


(43) 


r,- — 


2m M 


(44) 


That is, at a hnite separation distance, the massive particles attain their maximum speeds 
V = c. We note that since the massive particles always evolve to this state regardless of initial 
separation (he., various amounts of supposed “potential energy” in the initial state with no 
kinetic energy) energy cannot be conserved in this system. States with different energies all 
evolving into a single high-energy state requires sources or sinks in energy. However, the 
low-energy, non-relativistic approximation is quite useful for describing the dynamics at low 
energies. Unfortunately, unlike the Coulomb repulsion, there exist no initial conditions for 
which the relativistic limit is avoided. 

C. Relativistic corrections 


The discrete relations in Eqs. (15)-(16) may be formally interpreted as differential equa¬ 
tions by applying the convention 


An dv 

W+l W t — , 

An dn 

with a similar relation for tn —)■ t{n). One then obtains 

dv ^ 2mc 

* = = 

an 


(45) 


(46) 


Note that this corresponds to an acceleration given by the force in Eq. (13) with corrections 


which are hrst-order in /3 = -. Unlike the non-relativistic limit, this acceleration explicitly 
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drops to zero as f —)■ c. Furthermore the explicit appearance of v in the force indicates a 
non-conservative nature to this force. Employing the chain rule as for the non-relativistic 
limit, we obtain the following equation for /d(r), 


ndl^ m , . 


Eq. (47) is separable and admits the closed-form solution 

m ^ 

rp \ M g P 

Jo) 1-/5' 
This may be inverted to yield a formula for v = fic 

(ro/r)™/^ 


v{r) = c 


1 + lE 


where W{z) is the Lambert-W function,^ dehned as the principal value of 


(47) 


(48) 


(49) 


(50) 


While the solution clearly satishes n(ro) = 0 and 


lim v{r) = c [1 -I- 0] = c, (51) 

r^oo 


the time required for this to happen (rigorously, for |c — n| < 2mc/M) is quite large, and 
unfortunately for the theory, this does not appear to agree very well with the simulation 
or exact solution (see Fig. |^, breaking down even before the non-relativistic approximation 
breaks down. There is an equally curious situation that occurs in electromagnetism. The 
general solutions to Maxwell’s equations for known sources rely on fairly complex expressions 
involving evaluation of physical quantities at retarded times. However, by expanding these 
expressions the lowest-order term is the instantaneous Coulomb term. This appears to be 
a rather deep result also showing up in quantum electrodynamics and quantum gravity.®^ 
The rehned approximation in this section is only part of the required correction to the non- 
relativistic limit, and some potentially “fortuitous” cancellation between this modihcation 
and the rest of the terms being neglected is required to obtain a result more accurate than 
the NR approximation. An example of this sort of fortunate cancellation from classical 
physics may be observed by considering the electric held due to an arbitrary conhguration 
of currents and charges, given by one of Jehmenko’s equations,^ 


E(r,f) 


1 

dvreo 


cR PR 


dr', 


(52) 
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FIG. 4: Comparison of Eq. (49) to the simulation/exact solution and non-relativistic limit. Surpris¬ 
ingly, this “improved” approximation breaks down long before the NR limit ceases to accurately 


describe the physics. 


where p is charge density, J is current density, R= |r — r'|, and the retarded time is given 
hj tr = t — R/c. Following an exercise in a popular text on electrodynamics,^ one may 


consider constant currents for which J = 0 (he., the third term in Eq. (52) disappears). In 
this case, a miraculous cancellation occurs, yielding 


E(r,t) = 


dvren 


R^ 


dr' 


(53) 


where the correction to the instantaneous Coulomb potential and the second term in Eq. (52) 
cancel perfectly. That is, despite the explicit appearance of corrections of order (3 and 
evaluation of functions at tr instead of instantaneous time t, the held turns out to be 
the instantaneous Coulomb-like contribution. The “relativistic correction” in the particle- 
exchange model appears to be analogous to evaluating the held at retarded times without 
including the additional corrections, resulting in a less-accurate result at short times. 


IV. SUMMARY 

In this paper we thoroughly examined a simple model for classical interactions through 
the exchange of mediating particles in which momentum conservation is enforced for each 
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collision. As demonstrated in simulations and analytic reasoning, the resulting interactions 
yield an effectively conservative theory at low energies with a 1/r force. The conservative 
approximation breaks down at high energies, and regardless of initial separation, the massive 
particles both eventually reach the maximum speed allowed by the physical mechanism of 
energy transfer within the system. 

The classical particle exchange analogy of ice skaters throwing a ball back and forth has 
typically been used as an illustration in public outreach presentations and in teaching, from 
general education science courses to introductory and advanced physics courses. However, 
the analogy has value as a physical system for students to investigate quantitatively. The 
phenomenon can be used in various contexts including homework, an in-class activity, a 
computational physics exercise, or assessment. Furthermore, it can be used at both the 
introductory and advanced level in the undergraduate curriculum. 

In introductory physics, students learning computational modelin^^ can investigate 
the phenomenon numerically. Derivation of the change in speed of a massive particle, 
6v = 2mc/M, using Conservation of Momentum (Eq. is a straightforward exercise in 
introductory physics. Students can also explore and describe the position-time and velocity¬ 
time graphs. Because position and velocity change abruptly, introductory students have 
the opportunity to fit a smooth function to values that change discretely. Furthermore, 


teachers can use this system to assess understanding of potential energy functions (Eq. 21) 
and conservation of energy. Having already studied systems of particles interacting via the 
inverse-square law, students can practice applying a similar analysis to the 1/r force, possi¬ 
bly preparing them for similar forces that arise in an E&M course. Finally, as shown in this 
paper, teachers can also use the system as an application in a junior/senior level course in 
mechanic^ or mathematical physics where students are expected to explore the model in 
its limit using more advanced computational and analytical techniques. 
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